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Abstract 

We find equivalent hypergeometric- and difference-equation-based formulas, Q{k,a) = 
G\{a)G 2 {oL), for k = —1, 0,1,..., 9, for that (rational-valued) portion of the total separability 
probability for generalized two-qubit states endowed with random induced measure, for which 
the determinantal inequality \p^'^\ > \p\ holds. Here p denotes a 4 x 4 density matrix and 
its partial transpose, while a is a Dyson-index-like parameter with a = 1 for the standard (15- 
dimensional) convex set of two-qubit states. The dimension of the space in which these density 
matrices is embedded is 4 x (4 -|- fc). For the symmetric case of A: = 0, we obtain the previously 
reported Hilbert-Schmidt formulas, with (the two-re[al]bit case) Q(0, ^) = (the standard two- 

qubit case) Q(0,1) = and (the two-quater[nionic]bit case) Q(0,2) = The factors G^ia) 

can be written as the sum of weighted hypergeometric functions pTp_i, p > 7, all with argument 
H = (|)^- We find formulas for the upper and lower parameter sets of these functions and, then, 
equivalently express G^ia) in terms of first-order difference equations. The factors Gj^(a) are equal 
to (||)“~^ times ratios of products of six Pochhammer symbols involving the indicated parame¬ 
ters. Some remarkable a— and ^-specific invariant asymptotic properties (again, involving || and 
related quantities) of separability probability formulas emerge. 

PACS numbers: Valid PACS 03.67.Mn, 02.30.Zz, 02.50.Cw, 02.40.Ft, 03.65.-w 
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I. INTRODUCTION 


In a previous paper pQ, a family of (a-specific) formulas was obtained for the (total) 
separability probabilities of generalized two-qubit states. Here, we examine a related quan¬ 
tity informing us of that portion of the separability probabilities associated with the de- 
terminantal inequality > \p\- Here, p denotes a 4 x 4 density matrix and its 

partial transpose, with a serving as a Dyson-index-like parameter. Of course, by the Peres- 
Horodecki conditions EE], a necessary and sufficient condition for separability in this 4x4 
case is that > 0, while |p| > 0 itself certainly holds, independently of any separability 
considerations. So, the total separability probability can clearly be expressed as the sum of 
that part for which \p^'^\ > |p| and that for which \p\ > \p^'^\ > 0. The former part will be 
the one of immediate concern here. 

To obtain the new formulas to be reported, we employ the Legendre-polynomial-based 
density approximation (Mathematica-implemented) algorithm of Provost |1], utilizing the 
previously-obtained moment formula [SI sec. H] (cf. 0) 



|p|)”)/(|p|‘) = 

X 


(q + |)„ (p + + 2 + 5 q)„ 

2‘» (fc + 3a + |)„ (2k + 6a + 

^ 1 k \ a.^ k \ 2 , 0 . 

\1 — n — — n — a,n + 2k + 2 + ba' / 


while in [Tj, the moment formula [71 sec. X.D.6] 

n\ {a + 1 )^ {2a + 1 ), 


(-2k - 1 - 6a)„ (a)„ (a + 


2«”(3a+i)„(6a+|) 


2n 


5-^4 


— —77-, a-[-1, 2cr-h 1 ^ 

1—?7,,n-l-2-|-5a,l — n — a, n — a’ 


had been utilized for the density-approximation purposes there. (These [random-induced 
measure 0 ] moment formulas had been developed based on calculations solely for the two- 
rebit [a = and two-qubit [a = 1] cases. However, they do appear, as well, remarkably, 
to apply to the two-quater[nionic]bit [a = 2] case [U]. No explicit formal extension of the 
Peres-Horodecki partial-transposition conditions 0 0 to two-quaterbit systems seems to 
have been developed, however [cf. [Tnl - fl^ ].) 

In [1], a-specific formulas {a = 1,2,..., 13 and y |, |, |) as a function of k for the total 
{\p^'^\ > 0) separability probabilities had been derived. Here, contrastingly, we will find 
fc-specific formulas {k = —1,0,1,..., 9) as a function of a for the indicated one (|p^^| > |p|) 


2 





of their two component parts. We utilize an exceptionally large number (15,801) num¬ 
ber of moments in the routine of Provost |1], helping to reveal-to extraordinarily high 
accuracy-the rational values that the corresponding separability probabilities strongly ap¬ 
pear to assume. Sequences (a = 1, 2,..., 30,...) of such rational values, then, serve as input 
to the FindSequenceFunction command of Mathematica to obtain the initial set of /c-specihc 
(hypergeometric-based) formulas for Q{k,a), which we, then, further manipulate. 


II. COMMON FEATURES OF THE /c-SPECIFIC FORMULAS 


For each k = —1, 0,1,..., 9, the FindSequenceFunction command yields what we can 
consider as a large, rather cumbersome (several-page) formula, which we denote by Q{k,a). 
It, in fact, faithfully reproduces the inputted rational-valued (separability probability) se¬ 
quences. This fidelity is indicated by numerical calculations to apparently arbitrarily high 


accuracy (hundreds of digits). (The difference equation results below [sec. IV will provide 
a basis for our observation as to the rational-valuedness of the separability probabilities.) 

In Fig. we show plots of Q{k, a) over the range a E [1,10], for k = —1,..., 9. For fixed 
a, we have Q(/ci,a) > Q{k 2 ,a), if ki > ^ 2 . In Fig. we show a parallel plot, exhibiting 
linear-like behavior, for \ogQ{k,a). 


A. Distinguished jFq function with 2 as an upper parameter 

In each of the eleven fc-speci£c formulas obtained, there is a distinguished ^Fq function, 
with the (omnipresent) argument of || = (|)^ (cf. [12] [HI Ex. 8.6, p. 159]), with 2 as one 
of the seven upper parameters (cf. [IS]). The lower (bottom) six parameters conform to the 
simple linear rule 

{^ 1 , & 2 , & 3 , & 4 , & 5 , = ( 1 ) 

2k 23 2k 5 2k 27 2k 29 2k 31 

+ io’“ + y + 2-“+y+ io’“+y + io’“ + y + io’“+*^+® 

The six upper parameters (aside from the 2 mentioned) can be broken into one set of two 
(summing to an integer), incorporating consecutive fractions having 6’s in their denomina¬ 
tors, and one set of four (also summing to an integer), incorporating consecutive fractions 
having 5’s in their denominators. 
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FIG. 1: Plots of Q{k,a) over the range a G [1,10], for k = —1,...,9. For fixed a, we have 
Q{ki,a) > Q{k2,a), if ki > k2. 


Log(Q{k,ci:)) 



FIG. 2: Plots of log Q{k, a) over the range a G [1,10], for k = —1,..., 9. For fixed a, log Q(ki, a) > 
logQ(k 2 ,a), if ki> k 2 . 
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( 2 ) 


For the set of two, the smaller of the two entries abides by the rule 



^6q; + 4 


k 

3 


+ 2 


k + 1 
3 



1 


where the (integer-valued) floor function is employed, and the larger entry by 



^6q; -|- 2 


k 

3 


+ 4 


k + 1 
3 



( 3 ) 


For /c = 1, for illustrative purposes, application of these two rules yields {o; a -|- y}, 
and for k = 5, we have {a -1- y, a -|- (We have that ui +U 2 is an integer. The sequence 
of those integers is reproduced in A004523 [“Two even followed by one odd”] and A232007 
[“Maximal number of moves needed to reach every square by a knight from a hxed position 
on an n X n chessboard, or -1 if it is not possible to reach every square”] in the On-Line 
Encyclopedia of Integer Sequences [https://oeis.org/ol.html].) 

For the complementary set of four upper parameters, the entries in order of increasing 
magnitude are expressible as 


and 



( 4 ) 


Uq — a + - [ 2 


For k = 1, for illustrative purposes, application of these four rules yields 

{a a q; 4- y, q; and for k = 5, we have jo; a 4- y, a + ^, a ^ 


B. Distinguished function with 1 as an upper parameter 

Each /c-specihc formula Q{k,a) we have found also incorporates a second jFq function 
(again with argument ||, which is, to repeat, invariably the case throughout this paper), 
having all its thirteen parameters equalling 1 less those in the function just described. (A 
basic transformation exists [consulting the HYP manual of C. Krattenthaler, available at 
www.mat.univie.ac.at, allowing one to convert the thirteen [twelve a-dependent parameters, 
plus 1] of this jFq function [that is, add 1 to each of them] to those of the other jFq one first 
described.) 
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C. The remaining p-Fp_i functions, all with p > 7. 


Now, all the remaining m hypergeometric functions yielded by the FindSequenceFunction 
command for each of the fc-speci£c cases possess, to begin with, the same seven upper 
parameter (2 plus those indicated in (|^, (|^ and (|^) and the same six lower parameters 
Q, as in the first jFq function detailed above (sec. II A). Then, the seven upper parameters 
are supplemented by from one to m 2 ’s, and the six lower parameters supplemented by from 
1 to m I’s. 

From fc = — 1 to fc = 9, the eleven observed values of m are 


{m_i, mo, mi, m 2 , m 3 , m^, ms, me, my, ms, mg} = (3, 5,5, 6 , 6 , 7, 9, 8 ,10,10,10}. (5) 


D. Large a-free terms collapsing to 0 

We now point out a rather remarkable property of the formulas yielded by the 
FindSequenceFunction command. If we isolate those (often quite bulky) terms that do not 
involve any of the hypergeometric functions described above, we hnd (to hundreds of digits 
of accuracy) that they collapse to zero. These terms, typically, do contain hypergeometric 
functions similar in nature to those described above, but with the crucial difference that 
the Dyson-index-like parameter a does not occur among their upper and lower parameters. 
Thus, we are left with formulas Q{k,a) that are simply sums of m^ -|- 2 weighted pFp_i 
functions (of a), p = 7,... ,7 + m^. 

III. DECOMPOSITION OF Q{k,a) INTO THE PRODUCT G}(a)G|(a) 

The formulas we have obtained Q{k, a) can all be written-we have found-in the product 
form G}( 0)^2 (a). The G^ia) factor involves the summation of the hypergeometric functions 
pFp-i indicated above, each such function weighted by a polynomial in a, the degrees of the 
polynomials diminishing as p increases. Let us first analyze the other (hypergeometric-free) 
factor Gi{a), primarily involving ratios of products of gamma functions. 
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A. Hypergeometric-function-independent factor Gi(a) 


Some supplementary computations (involving an independent use of the FindSequence- 
Function command) indicated that this (hypergeometric-free) factor might be written quite 
concisely as 

i ('^2 )q,_i ('»3)q-_i ('»4)a_i ('»5 )q,_i ('»6)q;_i . x 

(b,l_, (b,)^_, (b,l_, (bel_, ’ ^ ^ 

where the Pochhammer symbol (rising factorial) is employed. Note that Gi(l) = 1. 


B. Hypergeometric-function-dependent factor Gl(a) 

1. Canonical form 

In Figs. |3]j^ we show a ’’canonical form” we have developed for the factors (cf. 

[13 Fig. 3]). 


IV. EQUIVALENT DIFFERENCE EQUATION FORMS 


It further appears that all the factors {k = —1,0,1,... ,9) (Figs. 1-4,...) can be 

equivalently written as functions that satisfy hrst-order difference (recurrence) equations of 
the form 

Poia) + Pi{a)G2{a) + p2{a)G2{l -|- a) = 0, (7) 


where the p’s are polynomials in a (Figs, 7-12). This was established by yet another appli¬ 
cation of the Mathematica FindSequenceFunction command. We generated-for each value 
of k under consideration-a sequence (a = 1,2,..., 85) of the rational values yielded by the 
hypergeometric-based formulas for ^^(a), to which the command was then applied. While 
we have limited ourselves in the last six hgures to displaying our results for k = —1, 0,1, 2, 3 
and 4, we do have the analogous set of results in terms of the hypergeometric functions 


for the additional instances, k = 5, 6, 7, 8 and 9, and presume that an equivalent set of 
difference-equation results is constructible (though substantial efforts with k = 5 have not 
to this point succeeded). The initial points G|(l) in the six difference equations shown are- 

in the indicated order-| Mif! members of the sequence 

f 188373 1096583 6050627 160298199 13988600951 1 = 1 Rip^p nrp Rip rP^nPpRvP 

I 785726 ’ 4342170 ’ 22951470 ’ 586426690 ’ 49611697974 / • uiicoc me uiic icnpceLivc 


are 
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in[6]:= TradltlonalForm[kmlnusl] 

3ut[6]//T raditional Form= 

1 


1517923 887480 

(a + 1) (2 O' + 3) (10 a + 9) (10 a + 11) (10 a + 13) (10 a + 17) {a (10 a+ 7) (925 + 615 „ + 134 ) + 54 ) 

1 1 2 3 4 5 9 11 13 3 17 27 

7F6I I, a + —, a -\ —, QfH—, QfH—, QfH—, an—;aH-, a + I, a -\ -, a H-, a -\—, a H-; — 

6 5 5 5 5 6 10 10 10 2 10 64 


a + ■ 


a + ■ 




a + 


5U 


a + 


1 


6 A 22 487 761296 


-125 (5 a (25 a (296 a + 303) + 2258) + 


7 6 7 8 9 11 19 21 23 5 

7F61 2 , O' H—, a H—, a H—, a H—, aH—, a H-; a H-, a + 2 , a -t -, a H-, an —, a + 

6 5 5 5 5 6 10 10 10 2 


1 


22487761296 


7 6 7 8 9 

625 (75 a (148 a + 101) + 1129) sAl 2,2, a —, a H —, a —, a H—, a H—, 

6 5 5 5 5 


11 19 21 23 5 27 27 

a H-; 1, a H-, a + 2, a -, a H-, an—, a H-; — 

6 10 10 10 2 10 64 


1 


22487761296 


7 6 7 8 9 11 19 

15 625 (296 a + 101) qFsI 2,2,2, an— , a -1 — , a -< —, an—,aH—,aH-; 1, 1 ,q'H- , a + 2, 

6 5 5 5 5 6 10 


21 23 5 27 273 1 

a H-, a H-, an—, a h-; — H- 

10 10 2 10 64j 11243 880648 


7 6 

578 125 10 F 91 2, 2, 2, 2, O' + -, a + - 
6 5 


7 8 9 11 19 21 23 5 27 27 

ffH—,Q'H—,Q'H—,aH-; 1, 1, 1,Q'H-,a + 2,Q'H-, a H-, an—, a H-; — 

5 5 5 6 10 10 10 2 10 64 


FIG. 3: G^Ha) 

separability probabilities Q{k, 1) themselves. We would like to extend this sequence suffi¬ 
ciently, so that we might be able to establish an underlying rule for it. (However, since the 
sequence is increasing in value, the Legendre-polynomial density-approximation procedure 
converges more slowly as a increases, so our quest seems somewhat problematical, despite 
the large number of moments incorporated [cf. [U App. II]].) If in the difference equation 







in[ 22 ]:= TraditionalForm [ kzsro ] 

Out[22J//T raditionalForm= 

-- il (a + 2)(2a + 3) (10 a + 13) (10 a + 17) (10 a + 19) 

779779185 6254401 

(10 a + 21) (a (5 a (25 a {2 a (740 a + 3119) + 10 313) + 208 403) + 410 694) + 63 000) 

1 2 3 4 5 7 13 3 17 19 21 27^1 

iFe\ 1, a -\—, q;h—, an—, an—, an—, an —;a -, a —, a h-, a h -, a + 2 , a h-; — + 

5 5 5 5 6 6 10 2 10 10 10 64j 


1 


5 a + 1 


1 

16875 laH— llaH— llan— l|aH— Man— l|aH— 


2 (5 a (5 a (5 a (75 a (296 a + 1089) + 115 606) + 389952) + 619097) + 362 847) 

/ 6 7 8 9 11 13 23 5 27 29 31 27 

jFel 2, a -I —, a H—, a H—, a H—, a H-, a H-; a H-, an—, a H-, a H-, a + 3, a -l -; — 

1 5 5 5 5 6 6 10 2 10 10 10 64 

f 6 7 

5 (30 a (5 a (50 a (370 a + 1089)+ 57 803) + 129 984)+ 619097) sFy 2, 2,a + -a + 

I 5 5 

8 9 11 13 23 5 27 29 31 27 

a -I —, a H—, a H-, a H-; 1, a H-, an—, a H-, a H-, a + 3, a -l -; — 

5 5 6 6 10 2 10 10 10 64 

f 6 7 8 9 

100 (5 a (25 a (1480 a+ 3267)+ 57 803)+ 64992) 9 F 8 2,2, 2, «+ a + -, a + -, a + -, 

I 5 5 5 5 

11 13 23 5 27 29 31 27 '1 

a H-, a H-; 1, 1, a H-, an—, a H-, a H-, a + 3, a -1 -; — + 

6 6 10 2 10 10 10 64j 


6 7 

125 (150 a (740 a + 1089) + 57 803) loFgl 2, 2, 2, 2, a + -, a + -, a + 


8 9 11 

—, a H—, a H-. 

5 5 6 


13 23 5 27 29 31 27'| 

a H-; 1, 1, 1 ,q;h- , a -1 —, an- , a -1 - , a + 3, a -1 -; — + 

6 10 2 10 10 10 64j 

6 7 8 9 11 13 

3750 (1480 a + 1089) iiFiol2, 2, 2, 2, 2, an —, a -1 —, a -1 —, an —, a -1 - , a -1 -; 1, 1, 1, 1, 

5 5 5 5 6 6 


23 


27 


29 


31 27^ 


a H-, a H—, a H-, a H-, a + 3, a -1 -; — + 925000 12 F 11 2, 2, 2, 2, 2, 2, a H—, a H— 


10 


10 


10 


10 64 ; 


8 9 11 13 23 5 27 29 31 27 

a -I —, an —, a -1 - , a -1 -; 1, 1, 1, 1, 1,q:h- , a -1 — , a -1 - , a -1 - , a + 3, a -1 -; — 

5 5 6 6 10 2 10 10 10 64 


Fig. 4 ; 62° (a) (the Hilbert - Schmidt case) 
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in[23]:= TradltlonalForm [ kplus 1 ] 

Out[23J//T raditionalForm= 

-- il (a + 3) (2 a + 5) (10 a + 17) (10 a + 19) (10 a + 21) (10 a + 23) 

4288 785 520 939 9201 

{a {a {a (20 a (5 a (740 a + 5783) + 91 541) + 3013 197) + 2724024) + 1 284280) + 246960) 

4 5 7 6 7 8 17 19 21 23 5 27^1 

iFe\ l,aH—, an—, an—, an—, an—, an—;aH-, a h-, a h-, a h -, an—, a + 3; — + 

5 6 6 5 5 5 10 10 10 10 2 64j 

16 875 la + -1 fa+ -1 fa+ -1 (a+ -1 fa+ -lfa+ - 


( 9 11 

(a (a (20 a (25 a (888 a+ 5783)+ 366 164)+ 9039591)+ 5448048) + 1284280) 7 F 6 2, a+ a+ —, 

I 5 6 

13 11 12 13 27 29 31 33 7 27l 

a H-, a H-, a H-, a H-; a H-, a H-, a H-, a H-, an—, a + 4; — + 

6 5 5 5 10 10 10 10 2 64j 

( 9 11 

(a (40 a (25 a (1110 a + 5783) + 274623) + 9039 591) + 2724024) gF? 2, 2, a + -, a + —, 

I 5 6 

13 11 12 13 27 29 31 33 7 27l 

a H-, a H-, a H-, a H-; 1, a H-, a H-, a H-, a H-, an—, a + 4; — + 

6 5 5 5 10 10 10 10 2 64j 


(40 a (25 a (1480 a+ 5783)+ 183082)+ 3013 197) 9 ^ 8 ( 2 , 2, 2, a + a + —, a + —, 

I 5 6 6 

11 12 13 27 29 31 33 7 27l 

an-,aH-,aH-; 1, l,aH-,aH-,aH-, an-,aH—,a + 4; — + 

5 5 5 10 10 10 10 2 64j 

/ 9 11 13 11 12 

20 (25 a (2220 a + 5783) + 91 541) 10 F 9 2, 2, 2, 2, a + -, a + —, a + —, a + —, a + —, 

I 5 6 6 5 5 

13 27 29 31 33 7 27'| 

a H-; 1, 1, l,aH-,aH-,aH-, an-,aH—,a + 4; — + 

5 10 10 10 10 2 64) 

f 9 11 13 11 12 13 

100 (4440 a + 5783) iiFio 2,2,2,2,2,an—,aH-,aH-,aH-,aH-,aH-; 1, 1, 1, 

5 6 6 5 5 5 


11 


13 


I 


27 


29 


31 


33 


27) 


1, a H-, a H-, an-,aH-,aH—, a + 4; — + 74000 12F11 2,2,2,2,2,2,an—,aH-, 


10 


10 


10 


10 


64) 


11 


13 11 12 13 27 29 31 33 7 27 

an-,aH-,aH-,aH-; 1, 1, 1, 1, l,aH-,aH-, an-,aH-,aH—, a + 4; — 

6 5 5 5 10 10 10 10 2 64 


Fig. 5 : 62^ (a) 
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in[25]:= TradltlonalForm [ kplus2 ] 

Out[25J//T raditionalForm= 

-'- (2 (a + 4)(2a + 5) (10 a + 21) (10 a + 23) (10 a + 27) (10 a + 29) 

625122980476 3944001 

(a (a (a (5 a (10 a (40 a (370 a + 4501) + 911 539) + 25032907) + 202090226) + 192 332 891) + 100092606) + 
22004136) 

f 7 6 7 8 9 11 21 23 5 27 29 27 \ 

I 6 5 5 5 5 6 10 10 2 10 10 64j 


2 (a (a (5 a (25 a (112 a (185 a + 1929) + 911 539) + 50065 814) + 303 135 339) + 192 332 891) + 50046 303) 

( 13 11 12 13 14 17 

7F6 2 , a H-, a H-, a H-, a H-, a H-, a H-; 

1 6 5 5 5 5 6 

31 33 7 37 39 27 '1 

an- ,q;h- ,a-\ —,q:h- , a -,q: + 5; — + 

10 10 2 10 10 64j 

13 

(2 a (5 a (50 a (420 a (74 a+ 643)+ 911539)+ 75 098 721)+ 303 135 339)+ 192 332 891)8^712, 2, a + —, 

6 


11 12 13 14 17 31 33 7 37 39 27'i 

an- , a -,q:h- , a -, a H-; 1, a -,q:h- , a —, a -,q:h-,q; + 5; — + 

5 5 5 5 6 10 10 2 10 10 64) 

( 13 

2 (10 a (25 a (280 a (185 a + 1286) + 911 539) + 25032907) + 101045 113) gFg 2, 2, 2, a + —, 

I 6 

11 12 13 14 17 31 33 7 37 39 27\ 

a -,q;h -,aH- ,q;h-, a H-; 1, 1, a H-, a H-, an—, a H-, an-,a + 5; — + 

5 5 5 5 6 10 10 2 10 10 64j 

( 13 11 12 

5 (50 a (280 a (370 a + 1929)+ 911 539)+ 25032907) 10 F 9 2, 2, 2, 2, a + —, a + —, a + —, 

I 6 5 5 

13 14 17 31 33 7 37 39 27\ 

a -,q;h -,aH-; 1, 1, 1, a H-, a -,aH—, q;h -, an- ,q; + 5; — + 

5 5 6 10 10 2 10 10 64j 

( 13 11 12 13 

50 (1680 a (185 a + 643) + 911 539) nFio 2,2,2,2,2,an-,aH- , a -{ -,aH-, 

I 6 5 5 5 

14 17 31 33 7 37 39 27'| 

a -, a H-; 1, 1, 1, 1 ,q:h-,q:h-,q:h—, an- ,q:h-,q: + 5; — + 

5 6 10 10 2 10 10 64j 

( 13 11 12 13 14 

14000 (370 a + 643) 12 F 11 2, 2, 2, 2, 2, 2, a H-, a H-, a H-, a H-, a H , 

I 6 5 5 5 5 

17 31 33 7 37 39 27', 

a H-; 1, 1, 1, 1, l,aH- ,q;h -,aH—,aH- ,q;h -,a + 5; — + 

6 10 10 2 10 10 64j 

( 13 11 12 13 14 17 

740000 13 F 12 2, 2, 2, 2, 2, 2, 2, a H-, a H-, a H-, a H-, a H-, a H-; 1, 

I 6 5 5 5 5 6 

31 33 7 37 39 27 

1, 1, 1, 1, 1,q;h -,aH-,aH—, an-, an-, a + 5; — 

10 10 2 10 10 64 


Fig. 6 : 62 ^ (a) 
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DifferenceRoot [Function [-[Oa a}, 

{(1+a) (2+a) (3 + 2a) (5 + 2a) (9 + lOa) (ll + lOa) {13 + lOa) (17 + lOa) (19 + lOa) 
(21+ 10 a) (23+ 10 a) (27 + 10 a) (54 + 938 a + 5645 + 12 625 + 9250 a^) + 

(-1 517 923 887 480 (2 +a) (5 + 2a) (19 + lOa) (21 + 10 a) (23 + 10 a) (27 + 10 a)) 
Ga'^ (a) + (569 221 457 805 (1 + 5 a) (2 + 5 a) (3 + 5 a) (4 + 5 a) (1 + 6 a) (5 + 6 a)) 

Ga-" (1+a) == 0, Ga-" (1) == —}]] 

14 

Fig. 7 ; Ga“^ (a) 


DifferenceRoot[ 

Function[{Ga“, a}, {(2+a) (3+a) (3+ 2 a) (5 + 2 a) (13+ 10 a) (17+ 10 a) (19+ 10 a) 
(21+ 10 a) (23+ 10 a) (27 + 10 a) (29 + 10 a) (31 + 10 a) 

(63000 + a (410694 + 5 a (208 403 + 25 a (10313 + 2a (3119 + 740 a))) )) + 

(-389 889 592812720 (3 +a) (5 + 2 a) (23 + 10 a) (27 + 10 a) (29 + 10 a) (31 + lOa)) 
Ga® (a) + (146 208 597 304 770 (2 + 5 a) (3 + 5 a) (4 + 5 a) (6 + 5 a) (5 + 6 a) (7 + 6 a)) 

Ga° (1+a) == 0, Ga° (1) == —}]] 

33 

Fig. 8 ; Ga° (a) (the Hilbert - Schmidt casej 


DifferenceRoot[ 

Function[{GaS a}, { (3 +a) (4 +a) (5 + 2 a) (7+2 a) (17 + 10 a) (19 + 10 a) (21 + 10 a) 

(23+ 10 a) (27+ 10 a) (29 + 10 a) (31 + 10 a) (33 + 10 a) (246 960 + 

a (1284 280 +a (2 724 024 +a (3 013 197 + 20 a (91 541 + 5 a (5783 + 740 a)))))) + 
(-2 144 392760469 960 (4 + a) (7 + 2 a) (27 + 10 a) (29 + 10 a) (31 + 10 a) (33 + 10 a)) 
Ga^ (a) + 804 147 285 176 235 (4 + 5 a) (6 + 5 a) (7 + 5 a) 

45 

(8 + 5 a) (5 + 6 a) (7 + 6 a) Ga^ (1 + a) == 0, Ga^ (1) == -}11 

286 


Fig. 9 : Ga^ (a) 


DifferenceRoot[Function[{Sa^/ a}, 

{(4+a) (5+a) (5 + 2 a) (7 + 2 a) (21 + 10 a) (23 + 10 a) (27 + 10 a) (29 + 10 a) (31+ 10 a) 
(33 + lOa) (37 + lOa) (39 + 10 a) (22 004 136 + a (100 092 606 + a (192 332 891 + 

a (202 090 226 + 5 a (25 032 907 + 10 a (911 539 + 40 a (4501 + 370 a))))))) + 
(-312 561490238 197 200 (5 +a) (7 + 2 a) (31 + 10 a) (33 + 10 a) (37 + lOa) (39 + 10 a)) 
Ga^ (a) + (117 210 558 839 323 950 (6 + 5 a) (7 + 5 a) (8 + 5 a) 

a a 1553... 

(9 + 5 a) (7 + 6 a) (11 + 6 a)) Ga^ (1 + a) == 0, Ga^ (1) == -}] 1 

8398 


Fig. 10 : Ga^ (a) 
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DlfferenceRoot[Function[{Sa^z a}, 

{(5+a) ( 6 +a) (5 + 2 a) (7 + 2 a) (27 + 10 a) (29 + 10 a) (31 + 10 a) (33 + 10 a) (37 + 10 a) 

(39 + lOa) (41 + lOa) (43 + 10 a) (134 548 128 + a (471 120 306 + a (698 007 782 + 

a (566 336789 + 5 0 (54 233 749 + 50 a (305 531 + 8 a (5833 + 370 a)) )))) ) + 

(-1 926733 941 573222 600 (6 + 0 ) (7 + 2 o) (37 + 10 o) (39 + 10 o) (41 + 10o) (43 + 10o)) 
62 ^ (o) + (722 525 228 089 958 475 (8 + 5 o) (9 + 5 o) (11 + 5 o) 

a a 3073 

(12 + 5 o) (11 + 6 o) (13 + 6 o) ) Ga^ (1 +o) == 0, Ga^ (1) == -}] 1 

14 858 


Fig. 11 ; Ga^ (o) 


DifferenceRoot[ 

Function[{Ga^, o}, { (6 +o) (7 + o) (7 + 2 o) (9 + 2 o) (9 + 4 o) (29 + 10 o) (31 + 10 o) 

(33 + 10 0 ) (37 + 10 0 ) (39 + 10 0 ) (41 + 10 o) (43 + 10 o) 

(47 + 10 o) (175 452 420 + o (522 054 355 + o (656 629 192 + 

o (451 645 197 + 4 o (45 743 164 + 5 o (2 174 607 + 400 o (698 + 37 o))) )))) + 

(-41795 305 501 819 136 400 (7 + o) (9 + 2o) (39 + 10o) (41 + 10o) 

(43 + 10 o) (47 + 10 o)) Ga^ (o) + 

(15 673239 563 182 176 150 (11 + 5 0 ) (12 + 5 o) (13 + 5o) (14 + 5o) (ll + 6o) (13 + 6o)) 

A A 8348 ... 

G2* (1+0) == 0, Ga* (1) == -} 

37 145 

Fig. 12 ; Ga^ (o) 
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( 8 ) 


for k = —1 (Fig. 7), we replace G 2 ^(1) = by G 2 ^(1) = 0 , then we can add 

(i)„+. (ll)„+ir(a)r(a + 2) 
52055003r(5a)r (a + 1) T (a + |) 

to the a-specific values obtained from the so-modi£ed equation to recover the values gener¬ 
ated by the original k = —1 difference equation (Fig. 7). 

A. Polynomial coefficients in difference equations 

We have for five {k = —1,1, 2,3,4) of the six cases at hand (Figs. 7-12) the proportionality 
relation 

P2(a) oc - 1), (9) 

where the UiS (and his) are themselves functions of both k and a. The (symmetric/Hilbert- 
Schmidt) case k = 0 fails to conform to this relationship because a factor of (1 -|- 5a) is 
present in the right-hand-side of ([^, rather than (6-|-5a), as in the corresponding difference 
equation. Now, for all six displayed cases (including fc = 0), 

pj(a) cx nl,b,. (10) 

Further, for all six cases, the polynomial coefficients p^ia) are proportional to the product 
of a factor of the form 

(11) 

and an irreducible polynomial. These polynomials are, in the indicated order, 

9250a^ + 12625a3 + 5645a^ + 938a + 54, (12) 

185000a'^ + 779750a^ + 1289125a^ + 1042015a2 + 410694a + 63000, (13) 

74000a® + 578300a® + 1830820a^ + 3013197a® -1- 2724024a2 + 1284280a + 246960, (14) 

and (for k = 2) 

740000a^ + 9002000a® + 45576950a® + 125164535a^ -t- 202090226a® (15) 

+ 192332891a2 + 100092606a + 22004136. 

The irreducible polynomial for fc = 3 is also of degree 7, that is, 

740000a^ + 11666000a® + 76382750a® + 271168745a^ + 566336789a® (16) 
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+698007782a2 ^ 471120306 a + 134548128. 


For = 4, this auxiliary polynomial is now the product of (9 + 4a) times an irreducible 
polynomial of degree 7, that is, 

296000a^ + 5584000a® + 43492140a® + 182972656a^ + 451645197a® (17) 

+656629192a^ + 522054355a + 175452420. 

The coefficients of the highest powers of a in all six irreducible polynomials are factorable 
into the product of 37 and powers of 2 and 5. 


V. Prob(|pT^ > 0) ANALYSES 


Efforts of our to conduct parallel sets of (fc-specific) analyses to those reported above 
for total separability probabilities > 0), rather than for that component part of the 

probabilities satisfying the determinantal inequality \p^'^\ > \p\ have so far been unsuccess¬ 
ful, in the following sense. We have computed what appear to be appropriate sequences 
(a = 1,2,..., 74) of rational values for /c = 1 and (a = 1,2,..., 124) for k = 2, but the 
Mathematica FindSequenceFunction has not produced any underlying governing rules. (This 
can be contrasted with the results in [1], where such successes were reported in obtaining 
a-specific [\p^'^\ > 0] formulas [a = 1, 2,..., 13 and |, |, |].) 

In Fig.[^ we plot the logs of these k = 1 seventy-four total separability probabilities (based 
on a = 1,..., 74). A least-squares linear fit to these points is —0.878482a — 0.362781, while 


in Fig. 10, we show (based on a = 1,..., 124) the k = 2 counterpart, with an analogous 
fit of —0.871033a -|- 0.351201. Although the slopes of these two linear fits are quite close, 
the y-intercepts themselves are of different sign. The predicted probabilities at a = 1, the 
first of the fitted points, are 0.289019 and 0.602955, respectively. In statistical parlance, the 
” coefficients of determination” or for the two linear fits to the log-plots are both greater 
than 0.99995. Further, sampling at a = 1, 51,101,..., 1451, we obtained an estimated, 
again, very-well fitting line of —1.4754 — 0.86417a. 
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Log[prob(|p'''^| > 0)] 



FIG. 9: Plot of logs of separability probability (|p^^| > 0) for random induced measure with k = 1. 
A least-squares linear fit to these 74 points is —0.878482 q; — 0.362781. 

Log[prob(|pPT|>0)]}] 



FIG. 10: Plot of logs of separability probability (|p^^| > 0) for random induced measure with 
k = 2. A least-squares linear fit to these 124 points is —0.871033a -|- 0.351201. 
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A. Asymptotic properties 


1. k-specific prob{\p^'^\ > 0) formulas 

Dunkl, on the basis of our k = l,a = 1, 51,101,..., 1451 analysis just above (and its 
companions), did advance the bold and (certainly, in our overall analytical context) elegant 
hypothesis of a fc-invariant (a —)■ oo) slope equal to log || ~ —0.8630462173553, which does 
seem quite consistent with the numerical properties we have observed (that is, with the di¬ 
rection in which the estimates of the slope tend as the number of points sampled increase). 
As further support, we obtained for a /c = 2, a = 1,49, 73,..., 1465 analysis, a slope estimate 
of -0.864025, again converging in the direction of log ||. (Let us remark, regarding the gen¬ 
eralized two-qubit version of the [simpler, lower-dimensional] X-states model [6l[l6l|T7j, that 
C. Dunkl has been able to show that the slope of a [now, log-log] plot of log(prob(|p^^| > 0) 
vs. logo tends to —as a —)■ oo.) 


2. a-specific prob{\p^'^\ > 0) formulas 


These interesting observations led us to reexamine, for their asymptotic properties, the 
’’dual” a-speci£c formulas reported in pQ. We now hnd-through analytic means-that for 
each of a = 1,2, 3,4 and that as fc —)■ oo, the ratio of the logarithm of the 

{k + l)-st separability probability to the logarithm of the fc-th separability probability is p. 
(Presumably, the pattern continues for larger a, but the required computations have, so far, 
proved too challenging.) For example, for a = we have for the two-rebit total separability 
probability, as a function of k, the formula [H eq. (4)] 


prebit _ 


A’^+^Sk + 15)r(A; + 2)r {2k + |) 
0fr(3fc + 7) 


(18) 


In Fig. 11, we show a plot of log(—(logP^®^**)) vs. k. The slope of a least-squares-htted line 
based on the 200 points is -0.523280, while log || ~ —0.523248. (As we increase a from 
but hold the number of points constant at 200, the approximation of the slope to this value 
slowly weakens.) 
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Log[-Log[prob(|pP'r|>0)]] 



FIG. 11: Plot of log(—(logP^®^**)) vs. k. The slope of a least-squares-fitted line is -0.523280, while 
log ^ « -0.523248. 

3. k-specific prob{\p^'^\ > |p|) formulas 

Now, as concerns the eleven (/c = —1,0,1,... ,9) formnlas for probdp-^^l > |p|), which 
have been the principal focns of the paper, we have compnted the ratios of the probability 
for a = 101 to the probability for a = 100. These ranged from 0.419810 {k = —1) to 
0.4204296 (fc = 9). Let ns note here that || ~ 0.421875. 

f. a-specific proh{\p^'^\ > |p|) formulas 

We had available a = |, 1 and 2 computations for k = 1,... ,40 for this scenario. We 
found that, for each of the three values of a, we could construct strongly linear plots-with 
unit-like slopes between 1.00177 and 1.00297-by taking k times the ratio {R) of the {k + 1) 
separability probability to the fc-th separability probability. (From this, it appears, simply, 
that i? —>■ 1, as fc —)■ oo.) 

5. ’’Diagonal” a = k prob{\p^'^\ > |p|) formulas 

For values a = k = 1,..., 50, we were able to construct a strongly linear plot by- 
similarly to the immediate last analysis-taking k = a times the ratio of the {k + 1) = (a-f 1) 
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separability probability to the k = a-th separability probability. Now, however, rather than 
a slope very close to 1, we found a slope near to one-half, that is 0.486882. The {k = a = 0)- 
intercept of the estimated line was 0.894491. 


VI. ’’CONCISE FORMULAS” 


Let us also remind the reader of the interesting ’’concise” (Hilbert-Schmidt [k = 0]) 
generalized two-qubit result-applying Zeilberger’s (’’telescoping”) algorithm na -of Qing- 
Hu Hon, reported in [151 eqs. (l)-(3)]. This-in our present notation-takes the form (cf. 


Figs. 5, 9) 

Q(0,a) = S“o/(a4-0, 

(19) 

where 

r,,. .run , 9(a)2-‘‘“-®r(3a + |)r(5a + 2) 

/(a)=Q(0,a) Q(0,a + 1)= ..xp/o ^ ^ l3^ ’ 

6r(a -I- l)r(2a -|- 3)r(5a -|- y) 

(20) 

and 




q{a) = 185000a® + 779750a^ + 1289125a® + 1042015a^ + 410694a + 63000 = (21) 

a (^5a (25a (2a(740a + 3119) + 10313) + 208403^ + 410694^ + 63000. 

We divide the originally reported formula by one-half, since we have moved here from the 
{k = 0) Hilbert-Schmidt > 0 original scenario to its \p^'^\ > |p| counterpart. Using 
our earlier results above, Hon has been able to construct the k = 1 analogue of the ’’concise 
formula”, 

Q(l,a) = E~o/(a-f-0, (22) 


where 


/(«) = 


q{a) (27)“r(5a)r(a + |)r(a + |) 


(50000)“ r(a)r (a + ^) r (a + |) r (a + |) r (a + I) r(2a + 5) 


(23) 


and 


g(a) = 


Ovr 


-(5a + l)(5a + 2)(5a + 3)x 


1000000 

(74000a® + 578300a® + 1830820a^ + 3013197a® -7 2724024a2 + 1284280a + 246960) . 


(24) 
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